Viscous Flow and Jump Dynamics in Molecular Supercooled Liquids: 
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The rotational dynamics of a supercooled model liquid of rigid A-B dumbbells interacting via a 
L-J potential is investigated along one single isobar. The orientation correlation functions exhibit 
a two-step character which evidences molecular trapping. Trapping is also apparent in a suitable 
angular-velocity correlation functions. The rotational correlation n is found to scale over more than 
three orders of magnitude as (T — T c ) ' with 7 = 1.47 ±0.01 and T c the MCT critical temperature. 
Differently, 77 with I — 2 — 4 and the rotational diffusion coefficient D r manifest deviations. For 
0.7 < T < 2 good agreement with the diffusion model is found, i.e. n = 1/1(1 + l)D r . For lower 
temperatures the agreement becomes poorer and the results are also only partially accounted for by 
the jump-rotation model. The angular Van Hove function evidences that in this region a meaningful 
fraction of the sample reorientates by jumps of about 180° . The distribution of the waiting-times in 
the angular sites cuts exponentially at long times. At lower temperatures it decays at short times 
as with £ = 0.34 ± 0.04 at T = 0.5 in analogy with the translational case. The breakdown of 
the Debye-Stokes-Einstein is observed at lower temperatures where the rotational correlation times 
diverge more weakly than the viscosity. 

PACS numbers: 64.70.Pf, 02.70.Ns, 66.20.+d , 66.10.-x 



I. INTRODUCTION 

The relaxation and the transport properties of molec- 
ular liquids depend on both their translational and ro- 
tational motion. Since their mutual interplay cannot be 
neglected both dynamical aspects must be jointly consid- 
ered. If the liquids are supercooled or supercompressed 
the overwhelming difficulties to molecular rearrangement 
are expected to enhance the role of the rotational de- 
grees of freedom and, more particularly, the rotational- 
translational coupling as noticed by experiments 
theory and numerical J&j f^] studies and discussed 
in a recent topical meeting [|13fl . 

Molecular-dynamics numerical ( MD ) studies provided 
considerable insight in supercooled liquids during the last 
years p4| . However, the question of rotational dynam- 
ics seems to have been partially overlooked since most 
studies dealt with atomic systems where rotational dy- 
namics is missing. Notable exceptions addressed the is- 
sue in model systems of disordered dipolar lattice ||, 
biatomic molecules [pJp| , pT| and well studied glassformers, 
e.g. CKN |5j, 0TP @0 and methanol JI| . The 
case of supercooled water was investigated in detail [[l9| . 
Studies of plastic crystals and orientational glasses ( i.e. 
no allowed translation ) are also known . 

We have recently presented numerical results on the 
the translational motion of a supercooled molecular 
model liquid |22| ( hereafter referred to as I ). The present 
paper wishes to complement I by extending the analysis 
to rotational degrees of freedom. As in I one issue is the 
detection and characterization of jump dynamics. It is 
found that rotational jumps are fairly more frequent than 



translational ones in the present system. This makes it 
easier their study. The occurrence of jumps poses the 
question of the coupling of the molecular reorientation 
with the shear viscous flow. This is the second issue ad- 
dressed in the paper. Jump dynamics may take place in 
the absence of any shear flow. Nonetheless, shear motion 
may favour jumps over energy barriers pq] . The ques- 
tion is of relevance in that the experimental situation 
is rather controversial. For macroscopic bodies hydrody- 
namics predicts that the reorientation is strongly coupled 
to the viscosity 77 according to the Debye-Stokes-Einstein 
law ( DSE ), t,D~ x cx 77, where r and D r are the rota- 
tional correlation time and diffusion coefficient, respec- 
tively p^-p6]|. DSE is quite robust. In fact, the coupling 
of the reorientation to the viscosity is usually found even 
at a molecular level if the viscosity is smaller of about 
1 — 10 Poise. At higher values DSE overestimates the cor- 
relation times of tracers in supercooled liquids according 
to time- resolved fluorescence |27|,^8|| and Electron Spin 
Resonance ( ESR ) studies P,| 29| ,|30|. On the other hand, 
photobleaching 0] and NMR yf studies found only small 
deviations from DSE even close to T g . Interestingly, ac- 
cording to ESR studies in the region where tracer reori- 
entation decouples by the viscosity, an ESR study evi- 
denced that it occurs by jump motion [ pl|]32[ . 

The paper is organized as follows. In section [n] details 
are given on the model and the simulations. In Sec. |III 
and 



IV 



the results are discussed and the conclusions are 



summarized, respectively. 
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II. MODEL AND DETAILS OF SIMULATION 

The system under study is a model molecular liquid 
of rigid dumbbells fici|-|l2]| . The atoms A and B of each 
molecule have mass m and are spaced by d. Atoms on 
different molecules interact via the Lennard- Jones poten- 
tial: 

V a p{r) = 4e Q/3 (a a/3 /r) 12 - (cr Q(3 /r) 6 , a,f3 e {A,B} 

(1) 

The potential was cutoff and shifted at r cuto ff — 
2A9(Jaa- Henceforth, reduced units will be used. 
Lenghts are in units of o aa, energies in units of caa and 

/ ^2 N.l/2 

masses in units of m. The time unit is ( maAA ] corre- 

V *AA J 

sponding to about 2ps for the Argon atom. The pressure 
P, temperature T and shear viscosity rj are in units of 
ZAA/&AA i e AA/k B and sJmt A Al '&\a-> respectively. 
The model parameters in reduced units are: a aa — 

GAB = 1.0, B B = 0.95, 6AA = £AB = 1.0, t B B = 0.95, 

d = 0.5, niA = tub = m = 1.0. The <jaa and <jbb values 
were chosen to avoid crystallization. The sample has N — 
N a ti '2 = 1000 molecules which are accommodated in a 
cubic box with periodic boundary conditions. Further 
details on the simulations may be found in I. 

We examined the isobar at P — 1.5 by equilibrating 
the sample under isothermal-isobaric conditions and then 
collecting the data by a production run in microcanonical 
conditions. The temperatures we investigated are T = 
6, 5, 3, 2, 1.4, 1.1, 0.85, 0.70, 0.632, 0.588, 0.549, 0.52, 0.5. 



III. RESULTS AND DISCUSSION 

This section will discuss the results of the study. We 
characterize the correlation losses of the system by inves- 
tigating several rotational correlation functions. Then, 
the related correlation times and transport coefficients 
will be presented. The presence of rotational jumps will 
be evidenced and their waiting time distribution will be 
discussed. Finally, the decoupling of the transport and 
the relaxation from the viscous flow will be presented. 



A. Correlation Functions 

The rotational correlation loss is conveniently pre- 
sented by suitable correlation functions. We study the 
dynamics of both the orientation and the angular veloc- 
ity of the dumb-bell. 



1. Orientation 



The rotational correlation functions are defined: 



^(*) = ^E< i ^( u *( t )- u *(°))) 

i=l 



(2) 



Uj(i) is the unit vector parallel to the axis of the molecule 
i at time t and Pi(x) the Legendre polynomial of order 
It is worth noting that C\ and Ci are accessible to sev- 
eral experimental techniques, e.g. dielectric spectroscopy, 
NMR, ESR, light and neutron scattering. 
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FIG. 1. Correlation functions C\ and C2 plotted for all the 
temperatures investigated. 

Fig. H shows C\ and C2. At high temperature and 
short times damped oscillations are present. They are 
typical features of free rotators in gas-like systems |^0|,[33| . 
At lower temperatures and intemediate times C\ and C2 
exhibit a wide plateau which evidences the increased an- 
gular trapping. At longer times the decay is fairly well 
described by the stretched exponential ~ exp[— (t/r) 13 ] 
with t = 62.9, (3 = 0.70 for I = 1 and r = 96.1, (3 = 0.60 
for I = 2 at T = 0.5. Even if C\ and C2 are quite similar 
two differences must be noted. First the plateau is lower 
at I = 2 than at I = 1. Second, at lower temperatures 
C2 vanishes at longer times than C\. The first feature 
is understood by noting that the oscillatory character of 
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P/(u,(t) • Uj(0)) with respect to the angle between Uj(i) 
and Uj(0) increases with I. Then, by increasing I even 
random angular changes with small amplitude occurring 
at short times affect the decay of Q. The second fea- 
ture is due to the fact that, as it will be shown later, 
molecules undergo frequent 180° flips at lower temper- 
atures. Due to the nearly head-tail symmetry, the flips 
reverse the sign of P;(uj(t) • Uj(0)) if I is odd whereas no 
change takes place if I is even. Then, they mainly affect 
the decay of Ci with odd I values. 





t 

FIG. 2. Correlation functions C3 and C4 for all the tem- 
peratures investigated. 

Fig. H plots the functions C3 and C4. The discussion 
is similar to the case I = 1,2. We note that, as expected, 
the plateau decreasess by increasing I, On the basis of 
the above discussion C3 should vanish before C4. How- 
ever, on increasing I, the effect of the head-tail symmetry 
is partially masked by the increased oscillatory character 
of the Legendre polynomials, which yields a larger sen- 
sitivity to small-angle reorientations. Similarly to Ci,2 
at longer times the decay is fairly well described by the 
stretched exponential ~ exp [-{t/rf] with t = 32.4, 
= 0.60 for I = 3 and r = 27.9, f3 = 0.47 for I = 4 



at T = 0.5. We notice that the stretching parameter 
decreases with increasing I. 

Fig. |^ compares the four correlation functions Ci with 
I = 1 — 4 at T = 0.5 are shown. Both the larger corre- 
lation loss at short times at larger I values and the odd- 
even effect on the long-time decay of the correlations are 
evidenced. 



u 




FIG. 3. Correlation functions Ci,Ci, C3 and C4 at T = 0.5. 
Ci decay more at short times by increasing I. The curves 
with I = 1, 3 decay faster at long times than the curves with 
I = 2,4, respectively, due to the partial head-tail symmetry 
of the molecules. 



2. Angular velocity 



For a linear molecule the angular velocity is: 



w = uxu 

A set of correlation functions is defined as: 

JV 

A" 



1 N 

*l(*) = -^S<fl(«»^(*))>. l>h 



(3) 



(4) 



i=l 



cei(t) is the angle between w<(t) and 07,(0). In particular, 
for I = 1 one has: 



(M 2 > 



(5) 



which is the usual correlation function of the angular ve- 
locity. In fig. ^ ^1 and ^2 are drawn for all temperatures 
investigated. \&i decays fast and increasing T slows down 
the decay ( note the difference with the orientation case 
). In particular, in the free-rotator limit is a constant. 
At lower temperatures shows a negative part at short 
times which evidences a change of sign of uj. This must 
be ascribed to the collisions on the rigid cage trapping 
the molecule ( see figs, [l] and ^ ). An analogous effect was 



3 



also noted for the linear velocity correlation function in 
I. More insight on the rotational trapping may be gained 
by inspecting ^ in fig^. No significant differences be- 
tween ^2 and ^1 are seen at high temperature. At lower 
temperatures, after a similar ballistic initial decay of 
and ^2, the latter slows down when $2 ~ 0.25. The 
long-living tail which shows up is interpreted by noting 
that at lower temperatures, after the ballistic regime the 
angular velocity is approximately trapped in a circle ( see 
eqj3| ). An elementary calculation shows that vE'i(t) = 
whereas ^(i) = 0.25 as long as the trapping is effec- 
tive pQ] . When the molecular rearrangement allows the 
orientation relaxation, the angular velocity tends to be 
distributed over a sphere and ^2 vanishes approximately 
as C2. 




A- = 5 



uj(0) ■ ui(t))dt 



(0) 



From a computational point of view the evaluation of 
the above integral is delicate and it is more convenient 
to evaluate D r via the Einstein relation 



D r = lim — - 

t^oc At 

R r is the mean squared angular displacement : 

N 

X 



1 N 



where </>i(£) is 



<j>i{t) - 4>i(0) = A&(i) = / Ui(t')dt' 



(7) 



(8) 



(9) 



In Fig. || R r (t) is shown. The plots are qualitatively 
similar to the mean squared translational dislaccmcnt ( 
see I ). At short time the motion is ballistic. At inter- 
mediate times and lower temperatures a plateau shows 
up. It signals the increasing trapping of the molecular 
orientation due to the severe contraints on the structure 
relaxation. At longer times the reorientation is diffusive 
according to eq.[?]. By comparing R r (t) with the trans- 
lational mean square displacement ( see fig. 3 of I) it is 
seen that the angular trapping is weaker than the one af- 
fecting the center-of-mass motion since the subdiffusivc 
intermediate regime is less pronounced and extends less 
on the time scale. 
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FIG. 4. Correlation functions of the angular velocity ^1 
and $2 plotted for all the temperatures studied. 

B. Diffusion coefficient and relaxation times 

The rotational diffusion coefficient of a linear molecule 
may be defined by a suitable Green-Kubo formula in close 
analogy to the translational counterpart as |33] : 
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FIG. 5. Mean squared angular displacement for all the 
temperatures investigated. 



The rotational correlation times are defined as 



n = / Q(t)dt 

Jo 



(10) 
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Fig. H presents the T-dependence of 77, I = 1 — 4 and 
It is seen that a wide region exists where the above 
quantities exhibit approximately the same Arrhenius be- 
havior ( about 0.7 < T < 2 ). At lower temperatures the 
apparent activation energy of the rotational correlat ion 
times increase. In particular, as noted in Sec. Ill A 1 , t\ 
becomes shorter than T2 and a similar crossover is an- 
ticipated between T3 and T4 at temperatures just below 
0.5. Differently, the rotational diffusion coefficient D r 
exhibits the same activated behavior over a region which 
was shown to extend also below the critical temperature 
T c predicted by the mode-coupling theory ( MCT ) jio) . 
The decoupling of D r with respect to 77 may be antic- 
ipated by noting that the former is related to the area 
below ^i(t) ( eq.0 ) and the latter to the area below 
Ci(t) ( eq.O ). At lower temperatures ^i(t) vanishes 
faster so probing the fast dynamics of the supercooled 
liquid whereas the d ecay of Ci(t) slows down more and 
more ( see Sec. [II A ). Alternatively, it has to be noted 
that even in highly-constrained liquids small angular mo- 
tions which are unable to relax the orientation lead to a 
finite value of D r in view of eq. ^ pj. Such librational 
motions were detected in an MD study of OTP p6[ . 

Fig. | shows also the MCT analysis of the T- 
dependence of the correlation time and the rotational dif- 
fusion JUJ. According to MCT, both 77 and D r should 
scale as 



t^D- 1 cx (T — T c y 



(11) 



The underlying expectation on the scaling [IT] is that it 
should work with the same T c value for any transport 
coefficient and relaxation time. Differently, the physical 
meaning of T c could be weakened. In I it was shown 
that eq.|ll] fits the divergence of the translational diffusion 
coefficient D over four orders of magnitude with T c = 
0.458±0.002 and ^ D = 1.93±0.02 ( the data are partially 
shown in fig. |^ ) and that at lower temperatures the 
primary relaxation time r a cx D^ 1 . Fig. |^ shows that the 
scaling is also effective for t\ with 7 = 1.47±0.01 ( see 
also fig.0). However, meaningful deviations are apparent 
for / > 2 and D r . The deviations increase as I increases. 
Since the rotational correlation functions Ci are more 
and more sensitive to small-amplitude reorientations and 
then to small molecular displacements, the poorer scaling 
may be a consequence of the difficulties of the mode- 
coupling theories at short distances pl,pa|. 
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FIG. 6. Arrhenius plot ( top ) and MCT scaling analysis 
( bottom ) of the rotational correlation times 77, / = 1 — 4 
and the rotational diffusion coefficient D r . T c — 0.458. The 
dashed lines are guides for the eyes. The translational dif- 
fusion constant D is also drawn for comparison ( open dia- 
monds ). The continuous line is the best fit by using eq.ll 
with 713 = 1.93 ± 0.02. 
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FIG. 7. MCT scaling analysis of n. T c = 0.458. 

To characterize the reorientation process, we studied 
the quantity 1(1 + l)D r Ti and the ratio 1(1 + 1)ti/2t\. 
If the the molecule rotates by small angular jumps, or 
equivalently the waiting time in a single angular site is 
fairly shorter than the correlation time 77, the motion is 
said to be diffusive and both quantities are equal to 1 for 
any I value |^,^3) . The results are shown in fig.|[ Three 
regions may be broadly defined. For T > 2 the properties 
are gas-like and the rotational correlation times become 
fairly long. For 0.7 < T < 2 l(l+l)D r T t and l(l + l)n/2ri 
do not change appreciably and are in the range 1 — 2. For 
T < 0.7 the above quantities diverge abruptly. The rapid 
increase in the deeply supercooled regime demonstrates 
the failure of the diffusion model which in fact is expected 
to work only in liquids with moderate viscosity or if the 
reorientating molecule is quite large. If the assumption 
of small angular jumps is released and proper account of 
finite jumps with a single average waiting time in each 
angular site is made, the socalled jump-rotation model is 
derived [3S|. The main conclusion is that 77 is roughly 
independent of I. In fact, for I = 2 it is found that the 
quantity 1(1 + l)n/2n w 3.5 at T = 0.5 ( see fig.| ) 
and the jump-rotation model predicts a value of about 
3. However, at higher I values the comparison becomes 
much less favourable. For I — 3 1(1 + 1)ti/2t\ ~ 1.9 at 
T = 0.5 whereas the prediction is about 6. For I = 4 
1(1 + 1)77/2ti w 2.75 to be compared to the prediction 
is about 10. The failure of the usual simple rotational 
models is not unexpected. Their basic assumptions are 
rather questionable in supercooled liquids, e.g. the in- 
herent homogeneity of the liquid and the presence of a 
single time scale both leading to the simple exponential 
decay of the rotational correlation functions. 
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FIG. 8. Temperature dependence of the quantity 
1(1 + l)D r n for I = 1, 2, 3,4 (top) and the ratio /(/ + 1)t;/2ti 
for 1 — 2,3,4 (bottom). If the reorientation is diffusive both 
quantities must be equal to 1. If the reorientation is jump- like 



C. Jump rotation 

The inadequate description provided by the diffusion 
and the jump model calls for a further characterization 
of the rotational motion. To this aim we consider the 
self-part of the angular Van Hove function: 



iV 



G e s (9,t) 



AT sin0 



(12) 



Qi(t) is the angle between the molecular axis of 
the i-th molecule at the initial time and time t . 
l/2G e s (9, t)sin0d9 is the probability to have the axis of a 



G 



molecule at angle between 9 and 9 + d9 at time t with re- 
spect to the initial orientation. At long times G e s (9, t) = 1 
since all the orientations are equiprobable. 

In Fig. H the function G e s is plotted for different tem- 
peratures and several times. At higher temperatures , as 
the time goes by, the molecule explores more and more 
angular sites in a continuous way. Instead, at lower tem- 
peratures G e s exhibits a peak at 9 w 180° and interme- 
diate times signaling that the reorientation has a mean- 
ingful probability to occur by jumps. 



T = 0.588 





FIG. 9. The angular Van Hove functions at low tempera- 
tures. Top: T = 0.588, bold line t = 80. Middle: T = 0.520, 
bold line t = 200. Bottom: T = 0.5, bold line t = 280. 
Dashed line t — 1. The other continuous lines are plotted at 
intermediate times with equal spacing. 

The indications provided by the VanHove function con- 
cerning the presence of rotational jumps are confirmed by 
directly inspecting the single particles trajectories ( fig. 
|To| ). Similar findings were reported also in other stud- 
ies on dumbbells and CKN glassformers Jl0| , p^5[ . With 
respect to the translational counterparts, it must be 
pointed out that they are quite faster ( see I ) and more 
frequent ( about one order of magnitude). The higher 
number of rotational jumps is also anticipated by noting 
that, differently from the translational Van-Hove func- 
tion, the rotational one does exhibit explicit signatures 
of jump motion ( see I ). 
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FIG. 10. Selected single-particle orientations evidencing 
the flips occurring at low temperatures. 



To characterize the jumps we studied the distribution 
^rot{t) of the waiting-time, namely the residence time in 
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one angular site of the unit vector being parallel to the 
axis of the i-th molecule. A jump of the i-th molecule is 
detected at t if the angle between Uj(to) and Uj(to+At*) 
is larger than 100° with At* = 24. To prevent multiple 
countings of the same jump, the molecule which jumped 
at time t is forgotten for a lapse of time At*. To minimize 
possible contributions due to fast rattling motion, each 
angular displacement is averaged with the previous and 
the next ones being spaced typically by 6 — 8 time units, 
depending on the temperature. The jump search proce- 
dure was validated by inspecting several single-molecule 
trajectories. The above definition of rotational jump fits 
well their general features, i.e. they are rather fast and 
exhibit no meaningful distribution of both the amplitude 
and the time needed to complete a jump ( sec fig.[Io] ). 
It is worth noting that in I it was found that the time 
needed to complete the translational jumps exhibits a 
distribution. The absence of a similar distribution for 
the rotational jumps points to a larger freedom of the 
latter. 
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FIG. 11. The rotational waiting-time distribution ip ro t(t) 
at low temperatures. 

Fig.|ll] shows ip ro t(t) at different temperatures. At 
T = 0.632 is virtually exponential. At lower tempera- 
tures deviations become apparent which are analyzed for 
T = 0.5 in fig. [l^. In I it was noted that the transla- 
tional waiting-time distribution may be fitted nicely by 
the truncated power law: 



<Kt) = [r(0T*] V 



-l e -t/T 



< £ < 1 (13) 



The best fit provided by eq.O is compared to the fits 
by using the stretched ( ezp] — (t/r)^] ) and the usual 
exponential functions in fig. n3. The better agreement 
of eq.[l3] at short times may be appreciated by looking at 
the residuals. 
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FIG. 12. Analysis of the rotational waiting-time distribu- 
tion iprot at T — 0.5. Top: comparison of the best fits with 
the exponential ( t = 102 ± 3 ), the stretched exponential ( 
stretching f3 = 0.78 ±0.02 and r = 63 ±4 ) and the truncated 
power law eq.13 ( £ = 0.34 ± 0.04 and r = 125 ± 3 ). The 
insert is a magnification of the short-time region. Notice that 
the fit with the power law is virtually superimposed to iprot- 
Bottom: residuals of the fits in terms of the truncated power 
law and the stretched exponential. 

The fractal behavior of both the translational and the 
rotational waiting-time distributions is an indication that 
the molecular motion at short times exhibits intermittent 
behavior. The issue in the framework of glasses has been 
addressed by several authors |p6[-p0[. The exponent £ 
of eq. 13 has a simple interpretation. If a dot on the 
time axis marks a jump, the fractal dimension of the 
set of dots is £. For £ < 1, it follows ip(t) oc t^ 1 at 
short times |36|,[4lJ in agreement with our results. If £ = 
1, the distribution of dots is uniform and ip(t) recovers 
the exponential form. This is expected beyond a time 
scale r and the exponential decay of ip ro t (t) at long times 
signals the crossover to the usual Poisson regime. It must 
be noted that the translational waiting time distribution 
at the lowest temperature ( T = 0.5 ) shows a weak 
tendency to vanish faster than eqjl^. A similar feature 
is not observed in ip ro t{t). 



D. Breakdown of the Debye-Stokes-Einstein Law 

For large Brownian particles the reorientation in a liq- 
uid occurs via a series of small angular steps, i.e. it 
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is diffusive. Hydrodynamics predicts that the diffusion 
manifests a strong coupling to the viscosity rj which is 
accounted for by the Debye-Stokes-Einstein law ( DSE ). 
For biaxial ellipsoids it takes the form [E5fl 



Di 



kT 



x,y,z 



(14) 



D 



x .y,z are the principal values of the diffusion tensor, k is 
the Boltzmann constant. The coefficients fXi depend on 
the geometry and the boundary conditions ( BC ). For 
a sphere with stick BC /J- X ,y.z = 6v, v being the volume 
of the sphere. For an uniaxial ellipsoid one considers 



D|| = D z and Dj_ = D x = A,. The case of stick BC 
can be worked analtycally [B4 |l 25 | . For slip BC numerical 
results for D± are known |26(] ( note that in this case 
the fluid does not exert torques parallel to the symmetry 
axis ). EqJIJ is sometimes rewritten in an alternative 
form in terms of proper rotational correlation times, e.g. 
for uniaxial molecules the equality r; = 1/1(1 + 1)-Dj_ 
holds. The new form is more suitable for comparison 
with the experiments since they do not usually provide 
direct access to the rotational diffusion coefficients. 

Irrespective of the heavy hydrodynamic assumptions, 
DSE works nicely even at a molecular level if the viscos- 
ity is not high ( 77 < lPoise ). Deviations are observed 
at higher viscosities for tracers in supercooled liquids by 
time-resolved fluorescence |27|j2g|1 and Electron Spin Res- 
onance ( ESR ) studies |p|, [29[ 30[ . On the other hand, 
photobleaching Q and NMRj3|, studies found only small 
deviations from DSE even close to T g . In all the cases 
known, DSE is found to overestimate the rotational corre- 
lation times since on cooling their increase is less than the 
one being exhibited by the viscosity. In this decoupling 
region ESR evidenced that the tracer under investigation 
rotates by jump motion |31|. 

Shear motion facilitates molecular jumps over energy 
barriers p^j . On the other hand, guest molecules may 
jump in frozen hosts in the absence of viscous flow. Since 
a meaningful fraction of the molecule in the system re- 
orientates by finite angular steps with intermittent be- 
havior, not quite expected in a liquid, it is of interest to 
investigate to what extent the reorientation is coupled to 
the viscous shear flow. 

The results are shown in fig. O by plotting the the 
quantity rj/XkT with X = D^ 1 ,l{l + l)n with 1 = 1- A. 
According to DSE it should be constant. The viscos- 
ity data were taken from I. At high temperatures the 
quantity approaches the value expected for stick BC. For 
T > 5 a tendency of ij/XkT for X = D r , t\ to increase 
is noted. However, at such temperatures the system 
manifests gas-like features ( see figs.[l],|],[| ) . On cool- 
ing r)/XkT increases. For intermediate temperatures the 
liquid properties are well developed, the system is dif- 
fusive ( 1(1 + l)nD r w 1, see fig.|| ) and -q/XkT has a 
value close to the DSE expectation with slip BC. No- 
tably, D r rj/kT remains close to this value in the wide 
interval 2 < T < 6. At lower temperatures rj/XkT di- 
verges. The stronger deviations are exhibited by D r and 



ti, the weaker ones by Ti. r 3 and r 4 track the behav- 
ior of Ti and T4, respectively, being the pair t\£ being 
affected by t he jum p motion much more than the pair 



t"2,4 ( see sec. Ill A 1 



rj increases of a factor of about 400 
between T = 1.4 and 0.5. The corresponding changes 
of D r r]/kT, rj/nkT and T)/T 2 kT are 41, 11 and 3.6, re- 
spectively. The discussion support the conclusion that 
the correlators being affected by rotational jumps, e.g. 
C\^(t) 1 yield correlation times fairly more decoupled by 
the viscosity. 

If one compares the changes of r\jT2kT to the ones 
drawn by ESR and fluorescence experiments in the re- 
gion T/T c 1.1 — 1.5, a broad agreement is found p7|,p[. 
These experiments found even larger values of ripr^kT 
on approaching T g . Instead, photobleaching studies de- 
tected changes of less than order of magnitude by chang- 
ing 7/ over about 12 orders of magnitude which are ex- 
tremely smaller than the present ones M. 
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FIG. 13. Plots of the quantity rj/XkT with 

X — D^ 1 ,1(1 + 1)77 over the overall temperature range ( top ) 
and in the supercooled regime ( bottom ). According to DSE, 
the quantity is a constant whose value depends on the hy- 
drodynamic boundary conditions. The superimposed dashed 
lines are the DSE expectations for stick and slip boundary 
conditions. Notice that at T = 1.4 rj/XkT depends little on 
X, signaling diffusive behavior. 

Inspection of figM suggests a way to reconcile, at least 
partially, the ESR ||, the photobleaching |J and NMR 
H studies on the glass former o-terphenyl (OTP). Both 
photobleaching and NMR measure T2 in a direct way. 
Instead, in the glass transition region the ESR lineshape 
depends in principle on several 77 and a model is needed 
to relate them to each other and lead to T2 |n]. The 
model is adjusted by fitting the teoretical prediction with 
the highly structured ESR lineshape. Fig.[l3] shows that 
the decoupling, as expressed by 77/77 fcT increases with 
I. Since the weight of the 77 to set the ESR lineshape 
is roughly comparable, one may anticipate that T2 SR 
may be underestimated at some degree around the glass 
transition. On the other hand, it must be pointed out 
that the rotational decoupling is observed up to 1.4T g in 
OTP where ESR yields T2 in a model-independent way 
||. Furthermore, the decoupling is evidenced also by 
fluorescence experiments which provide T2 in a model- 
independent way 27 1 . 



IV. CONCLUSIONS 

The present paper investigated the rotational dynam- 
ics of a supercooled molecular system. The study ad- 
dressed several general features and focussed on the char- 
acterization of the jump dynamics and the degree of cou- 
pling with the viscosity. 

The ensemble consists of rigid A-B dumbbells interact- 
ing via a L-J potential. All the properties were studied 
along the isobar P — 1.5. The time orientation correla- 
tion functions Ci (t) exhibit at high temperatures gas- like 
features. In the supercooled regime, after a first initial 
decay, a plateau is observed which signals the trapping 
of the molecule due to the increased difficulty of the sur- 
roundings to rearrange themselves. The plateau level de- 
creases by increasing the / rank of the correlation func- 
tion due to the larger sensitivity to small-angle librations. 
The long-time decay of Ci(t) is fairly well decscribed by 
a stretched exponential. The stretching increases from 
I = 1 ( /3 = 0.7 ) to I = 4 ( /3 = 0.47 ) at T = 0.5. The 
influence on the decay time of Ci ( and then on 77 ) of the 
partial head-tail symmetry of the dumbbells was noted. 

A set of angular- velocity correlation functions \&; was 
defined. ^1 decays faster and faster on decreasing T. 
Differently, ^2 develops a long-lasting tail which vanishes 
on the time scale of C%. 

The temperature dependence of the rotational correla- 
tion times 77 with 1 = 2—4 and the rotational diffusion co- 



efficient D r manifest deviations by the power-law scaling 
in T — T c , being T c the MCT critical temperature. There 
were ascribed to the difficulties which mode-coupling the- 
ories meet at short length scales ^J,^5[ . Remarkably, the 
scaling works nicely for n over more than three orders 
of magnitude with an exponent 7 = 1.47 ± 0.01. This 
parallels the scaling which was noted for the translation 
diffusion coefficient and the primary relaxation time T a 
in I. For 0.7 < T < 2 the quantities 1(1 + l)D r n and 
1(1 + 1)t; /2ti do not change appreciably and are in the 
range 1 — 2 in good agreement with the diffusion model 
which predicts that both of them equal to 1 indepen- 
dently of the temperature. For T < 0.7 the above quan- 
tities increase abruptly. The increase of the quantity 
1(1 + 1)77/271 is reasonably accounted for by the jump- 
rotation model for I = 2 p2]| . For higher I values the 
agreement becomes quite poor. Being t± and 72 mea- 
sured by most experiments p|-p|, p8p2] | this finding may 
account for the attention that the jump model has at- 
tracted during the last years . 

The analysis of the angular Van Hove function evi- 
dences that in this region a meaningful fraction of the 
sample reorientates by jumps of about 180°. The flips are 
rather fast and exhibit no meaningful distribution of both 
the amplitude and the time needed to complete a jump. 
Differently, it was noted in I that translational jumps re- 
quire different times to be performed. The absence of a 
similar effect for the reorientations indicates a larger an- 
gular freedom. This is also apparent by the larger num- 
ber of rotational jumps which are detected with respect 
to the translational ones. It is worth noting that the ease 
to jump does not lead to trivial relaxation properties, as 
signaled by the stretched decay of 61,2.3,4- 

We characterized the distribution i/v ot of the waiting- 
times in the angular sites. It vanishes exponentially at 
long times whereas at lower temperatures it decays at 
short times as i 5 " 1 with £ = 0.34 ± 0.04 at T = 0.5. 
Interestingly, the translational waiting-time distribution 
exhibits the same behavior ( see I ). The exponent for 
the translational case is £ = 0.49. We ascribe the power- 
law to the intermittent features of the motion in glassy 
systems 

The intermittent jump reorientation is fairly different 
from the motion in a liquid. Then, a decoupling from 
the viscous flow and the subsequent breakdown of the 
Debye-Stokes-Einstein is anticipated. Our study con- 
firms the breakdown and shows that the quantity rj/XkT 
with X = D^ x ,l(l + 1)77 and I — 1 — 4 diverges below 
T = 1. In particular, the correlators being affected by ro- 
tational jumps, e.g. Ci s(i), yield correlation times fairly 
more decoupled by the viscosity. A rather similar decou- 
pling was found in I for the product Drj, D being the 
translational diffusion coefficient. 

The decoupling of the molecular reorientation by the 
viscosity could be also anticipated by the observed ease 
to perform rotational jumps. The reduced tendency to 
freeze of the rotational degrees of freedom was pointed 
out by MD M and theoretical M studies. The former in- 



10 



vestigating the residual rotational relaxation in a random 
lattice with quenched translations and the latter predict- 
ing a hyerarchy for the glassy freezing, i.e. the rotational 
dynamics can never freeze before the translational dy- 
namics. The decoupling of the rotational motion of guest 
molecules from the viscous flow has been experimentally 
seen by by time-resolved fluorescence |27],^8| and Elec- 
tron Spin Resonance [p|, p9|j3C| l while photobleaching and 
NMR studies reported small deviations from DSE even 
close to T g HH . ft is worth noting that the decoupling of 
the translational diffusion from the viscosity and related 
phenomena as the so-called rotation-translation paradox 
have been ascribed to a spatial distribution of mobility 
and relaxatio n prop erties, so called dynamical hetero- 
geneities Pfl^5|-^9| . Their role will be addressed in a 
forthcoming study. 
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